ON THE MOTION UNDER FOCAL ATTRACTION IN 
A ROTATING MEDIUM 



Jorge Sotomayor 

Institute de Matematica e Estatfstica, Universidade de Sao Paulo 
Rua do Matao 1010, Cidade Universitaria 
05.508-090, Sao Paulo, SP, Brasil 
sotp@ime.usp.br 



Abstract. New results are established here on the phase portraits and bifur- 
cations of the kinematic model in system JTJ, first presented by H.K. Wilson in 
[3], and by him attributed to L. Markus (unpublished). A new, self-sufficient, 
study which extends that of [3] and allows an essential conclusion for the 
applicability of the model is reported here. 

1. Introduction 

Consider the following family of planar differential equations depend- 
ing on three real parameters [uj, v , R), with uj > 0, v > 0, R > 0, 

x = —uy + v- 

(1) 

V = UJX — V 



The solutions, also called orbits, of system (OQ) describe the motion 
of certain entities (such as particles or micro-organisms) attracted by 
a focus F = (R, 0) (such as a light source or a magnetic pole) toward 
which they move with velocity v; the motion takes place in a medium 
(such as a fluid) which rotates with angular velocity uj around a fixed 
point, located at the origin O = (0, 0). 

The focal point F = (R,0), where system ([T|) is undefined, will be 
regarded as a singular point, outside of which it is analytic. A point 
at which the components of the system vanish will be referred to as an 
equilibrium point. 

According to Wilson [3J, p. 297, this is a model suggested by L. 
Markus (see [3], p. viii,) for the motion of phototropic platyhelminthes 
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— light-seeking flatworms — swimming in a liquid filling a shallow cir- 
cular recipient with section 

G R = {x 2 + y 2 <R}. 

No printed reference source containing Markus' suggestion is given in 

i- 

Definition 1. Denote by W S (Q) the basin of attraction of an equilib- 
rium or singular point Q of This is the set of points po such that 
y(^Po) Q as t —>■ m + (po), the right extreme of the maximal interval 
of<p(t,po). The basin of repulsion, W U {Q), is defined analogously for 
t — > m_(p ), referring to the left extreme of the maximal interval. 

Notice that for Q = F the approach to F by orbits of ([T]) happens 
in finite time, that is, m + (po) or m_(po) is finite. 

In [3J, p. 298, using the Poincare - Bendixson Theorem and the 
Bendixson Negative Criterium, it is proved the existence and unique- 
ness of an equilibrium point P = Pf UjVi R) attracting all positively com- 
plete semi-orbits <p(t,p ) of system ([!]) starting at points p £ Gr \ F. 
Complete means that ip(t,p ) is defined for all t G [0, oo). Denote such 
set of points p by Gr i+ . 

For the terms in ODE's not defined here, besides [3J, the reader can 
profit from consulting Chicone jl] , among other more recent up to date 
books. 

Under the assumption to > v/R, P is 

(2) P {u>v , m = ((v/uj) 2 /R, yjR?-(v/u) 2 (v/u)R)). 

This corrects a misprint in [3J, p. 298. 

The theorem below improves the results on this subject outlined in 
[3J . An elementary proof will be given in section [2J 

Theorem 1. For alluo > 0, the region Gr\F is positively invariant, in 
fact the radial component of system (Q]j is negative on the complement 
of the closed disk C of center (R/2, 0) and radius R/2. 
It holds that 

1. For < co> < v/R, F is a global attractor: W S (F) =R 2 \F. 

2. For uj > v/R there is a unique hyperbolic attracting equilibrium 
P located at whose basin of attraction contains W U (F), the 
basin of repulsion of F, itself a regular analytic curve contained 
inG R \F. 

3. Also, W S (F), the basin of attraction of F is a regular analytic 
curve disjoint from Gr . 

In particular it holds that for u > v/R, 

g r , + = G R nw s (P) = G R \F 
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2. Proof of Theorem Q] 

Performing the change of variables and parameter rescaling 

(3) x = Rx, y = Ry, t = iR/v, u = Cuv/R 

and then removing the bars, obtain system ([T]) with R = v = 1 : 

,., 1 — x . y 

(A) x = -ujy + = , y = ujx =. 

Writing this equation in polar coordinates centered at F, given by 

(5) x = 1 — r cos 8, y = r sin 9, 

and multiplying both components by r, which amounts to rescaling 
again the time, obtain 

(6) 9 = — uj{r — cos 9), r = r(uj sin9 — 1). 

Derivation of z = x 2 + y 2 (which, in polar coordinates, is z = 1 + r 2 — 
2r cos 9) in the direction of equation (j4]) (which, in polar coordinates, 
has the same direction field as gives z' = — 2r(r — cos#). 

Clearly z' is negative above the curve r = cos 9, which is the polar 
expression for the border of the disk C. 

Calculation of the equilibrium of (jUJ) on the plane r > gives the 
point P whose polar coordinates (9p,rp) satisfy 

sin 9p = 1/uj, rp = cos9p = \Jl — {l/uj) 2 . 
The expression, (xp,yp), of P in cartesian coordinates follows: 

x P = l- y/l - 0M 2 fi - (1/uj) 2 = (1/u) 2 , y P = sjl - (l/wJVw, 

which in the original coordinates, before the change in (j3J), reproduces 
the expression in (|2"1). 

Calculation of the divergence, a, and Jacobian, 5, of (jSJ) gives 

a = — 1, 5 = —uj 2 + u) sin 9 + u) 2 r cos 9 + u) 2 cos 2 9. 

Evaluation at P, i.e. at (9p,rp), gives 5 — u 2 — 1. This shows that 
P is an attracting hyperbolic equilibrium point of node or focus type. 

Calculation of the discriminant, a 2 — AS, for values of oj > 1 gives 
that the transition of P from a node to a focus occurs at y/5/2, that is 
at u = \/o~v/2R, in the original parameters. 

This proves the first general assertion and item 1 in the theorem. 



4 



JORGE SOTOMAYOR 



For u > 0, system ([6]) has two equilibria; one is a hyperbolic sad- 
dle at S- = (— 7r/2,0) with unstable separatrix along the #-axis and 
attracting eigenvector parallel to (1, (2u> + 

The other equilibrium, located at S + = (vr/2, 0), is a hyperbolic node 
for uj < 1. For uj > 1, S + = (vr/2, 0) is a hyperbolic saddle with stable 
separatrix along the #-axis and with repelling eigenvector, for positive 
eigenvalue uj — 1, parallel to (1, (— 2u + 1)/uj). 

For the equilibrium at (vr/2, 0), the transition at uj = 1, for increasing 
co>, gives a cubic pitchfork node to saddle bifurcation. For > 1, 
bifurcate two attracting nodal equilibrium points, one to r > and the 
other to r < 0, which capture the unstable separatrices of the saddle. 
Only the first one, located at P, is seen in the original (x, y)-plane. 




Therefore, the slope of the unstable separatrix, at S + , is —2 + l/u 
and that of the stable separatrix at S_ is 2 + 1/uj. 

Comparison at S- of the slope the stable separatrix with the slopes 
of the borders of Gr (for R = 1), given by r = 2 cos#, which is 2, and 
of C, given by r = cos#, which is 1, leads to the location of W S (F) 
outside of Gr. 

The location of W U (F) between Gr and C follows from similar com- 
parison of slopes, taking into account that at S + the slopes of the 
borders of Gr (for R — 1) is —2, and of C is —1. 

See Fig. [fl in polar coordinates, and Fig. [21 in the original coordi- 
nates. 

By continuity, for any u > 1 (that is greater than v/R, in the original 
coordinates), the basin of repulsion of F is always contained in W S {P). 
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Focus (R, 0) 



Figure 2. Borders of G and C, Equilibrium Point P 
and Separatrix Directions at Focus F 

This, for uo near 1, is a consequence of the structure of the local saddle- 
node bifurcation. For large u this follows from Poincare-Bendixson 
Theorem and the Bendixson Negative Criterium [3]. 

From the fact that z' is negative outside the disk C follows that 
W S (P) coincides with R 2 \ W S (F). 

The two main conclusions in items 2 and 3 are established. 



3. Concluding Comments 

Item 3 is essential for the purpose of the biological model in [3j. In 
fact, from the knowledge of the correct location, at (J2]), of the equi- 
librium attracting all entities moving with velocity v, their separation 
from similar ones also contained in the positive invariant set Gr but 
having different characteristic velocities and, therefore, clustering at 
different points of Gr, for large t. From this information the removal 
of the entities for study in isolation could be implemented practically. 
However, no laboratory experiment report where the model has actu- 
ally been used is known by the author. 

Both Wilson [3] and Sotomayor [2], p. 264, the only written sources 
for system (pQ), overlooked the property in item 3 of the theorem. Its 
consideration was proposed to the author toward 1995 by Dan Henry 
(1944 - 2002). 
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